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摘 要 :基于 非 局 部 理论 和 分 数 阶 导 数理 论 , 研 究 上 履 夭 弹性 场地 土 的 地 震 放大 效应 。 利 用 Erin- 
gen 非 局 部 理论 考虑 土 体 颗粒 尺度 等 非 局 部 效应 的 影响 ,通过 分 数 阶 黏 弹 性 本 构 模 型 刻画 场地 土 的 
应 力 应 变 本 构 关系 ,建立 基于 非 局 部 理论 的 分 数 阶 黏 弹 性 场地 土 的 振动 微分 方程 ;考虑 分 数 阶 导 数 
各 的 性 质 和 黏 弹性 场地 土 的 边界 条 件 , 得 到 了 简 谐 地 震波 作用 下 黏 弹性 场地 土 的 位 移 和 剪 切 应 力 的 
解析 解 ,并 在 频率 域内 给 出 了 位 移 放 大 系数 和 应 力 放 大 系数 的 表达 式 ; 最 后 通过 数值 算 例 分 析 了 非 
局 部 效应 、 分 数 阶 导数 的 阶 数 和 土 体 黏 性 参数 等 对 黏 弹性 场地 地 震 放 大 效应 的 影响 。 数 值 分 析 结 
果 表 明 , 在 低频 时 位 移 放大 系数 和 应 力 放大 系数 随 频 率 变 化 曲线 存在 波动 ,高 频 时 逐渐 趋 于 稳定 ; 
非 局 部 效应 对 场地 土 位 移 放大 系数 的 影响 与 频率 有 关 , 对 应 力 放 大 系数 的 影响 较 大 ,在 研究 场地 土 
振动 效应 时 有 必要 考虑 土 体 非 局 部 效应 的 影响 ;分 数 阶 导数 的 阶 数 越 小 ,位 移 放 大 系数 和 应 力 放 大 
系数 随 频 率 变化 曲线 波动 越 大 ;场地 土 的 力学 性 质 对 场地 土 的 振动 效应 的 影响 较 大 ;上 履 场 地 土 的 
时 黏 性 对 位 移 放大 系数 的 影响 与 频率 有 关 , 高 频 时 , 土 体 莫 性 越 大 ,位 移 放 大 系数 越 大 ; 越 接近 基 洗 ， 
土 体 的 应 力 放大 系数 越 大 , 且 土 体 深度 对 应 力 放大 系数 的 影响 越 大 。 
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Abstract :Based on nonlocal theory and fractional derivative theory ,the seismic amplification effect of o- 
verlying viscoelastic soil is studied. The dynamic differential equations of fractional viscoelastic ground soil 
are established based on the nonlocal theory by considering the influence of nonlocal effects such as parti- 


cle size ,and the nonlocal effects is considered by using Eringen’s nonlocal theory , and the stress-strain 


收 稿 日 期 :2021-01-25 修 回 日 期 :2022-06-03 
基金 项 目 : 国 家 留学 基金 资助 项 目 (No.201908410098 ) ;河南 省 高 等 学 校 重 点 科研 资助 项 目 (No.21A560017 ) 


通信 作者 : 谢 顺 利 ,讲师 。E-mail:53845664@ qq. com 


引用 格式 : 谢 顺利 , 张 春 丽 , 周 阿 祥 . 基于 非 局 部 理论 的 分 数 阶 颖 弹性 场地 地 震 放 大 效应 分 析 [ 中 .应 用 力学 学 报 ,2022 ,39(4) : 726-731 ,741. 


XIE Shunli, ZHANG Chunli,ZHOU Axiang. Analysis of vertical dynamic amplification effect of fractional viscoelastic ground soil based on non- 


local theory[ J ]. Chinese journal of applied mechanics ,2022 ,39(4) :726-731 ,741. 


constitutive relation of ground soil is described by using the fractional derivative viscoelastic constitutive 
model. Considering the properties of fractional derivative and boundary conditions of viscoelastic ground 
soil] ,the analytical solutions of displacement and shear stress of viscoelastic ground soil under harmonic 
seismic wave are obtained ,and the expressions of displacement amplification coefficient and stress amplifi- 
cation coefficient are given in frequency domain. The numerical results show that the displacement amplifi- 
cation coefficient and stress amplification coefficient fluctuate with frequency at low frequency ,and tend to 
be stable at high frequency. The influence of nonlocal effect on the displacement amplification coefficient 
of ground soil is related to frequency ,and the influence on the stress amplification coefficient is greater. It 
js necessary to consider the influence of nonlocal effect when studying the dynamic effect of ground soil. 
The smaller the order of fractional derivative ,the more severe the fluctuation of curves of vertical displace- 
ment amplification coefficient and stress amplification coefficient varying with frequency. The mechanical 
properties of ground soil have a greater impact on the vertical dynamic effect of ground soil. The influence 


of the viscosity of the overlying soil on the displacement magnification coefficient is related to the frequen- 


cy, When the high frequency and the larger the viscosity of the soil ,the larger the displacement magnifica- 


= tion coefficient. The closer to the bedrock ,the greater the stress amplification coefficient of soil, and the 


CN larger the stress amplification coefficient changes with depth. 
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@ 〇 在 地 震波 作用 下 ,场地 土 对 地 震波 有 放大 和 滤 
波 作用 ,场地 条 件 和 场地 土 的 力学 性 质 对 场地 土 的 
地 震 放 大 效应 的 影响 重大 ,研究 不 同 场地 条 件 和 考 
虑 场地 土 各 部 分 区 域 对 地 震 放大 效应 的 影响 对 地 下 
结构 或 建筑 物 的 动力 响应 和 稳定 问题 具有 十 分 重要 
的 工程 应 用 价值 , 事 关 工 程 安全 问题 。Idriss 等 " 基 
玉 汪 维 剪 切 梁 模 型 研究 了 剪 切 模 量 为 常数 和 沿 深度 
授 乱 函数 变化 的 场地 土 的 水 平 振动 特性 和 地 震 反 
应 SBoer 等 2 利用 多 孔 介 质 理论 得 到 了 不 可 压 流 体 
多 卉 半空 间 一 维 动力 响应 的 解析 解 ,并 研究 了 平面 
波 的 传播 和 反射 间 题 ;Davis™“ 借助 等 效 剪 切 模 量 对 
剪 切 模 量 沿 深度 按 宕 函数 变化 的 场地 土 的 振动 问题 
进行 了 研究 ;高 玉 峰 等 5” 在 时 间 域 内 研究 了 成 层 地 
基 一 维 土 层 地 震 反应 问题 ,得 到 了 问题 的 解析 解 。 
栾 茂 田 等 研究 了 水 平成 层 匀 质 场地 自 振 特性 ,并 
给 出 了 时 域内 地 震动 力 响应 的 解析 表达 式 。 尚 守 平 
等 "在 一 维 波动 模型 的 基础 上 ,研究 了 坚 直 向 上 传 
播 的 稳 态 剪 切 地 震波 作用 下 剪 切 模 量 沿 深度 按 指数 
规律 变化 的 场地 土 的 横向 自由 振动 问题 ; 杨 吴 等 
在 多 孔 介质 理论 的 基础 上 ,借助 Crump 数值 反 演 研 
究 了 笑 弹 性 饱和 不 可 压 多 孔 介 质 层 的 一 维 动力 响应 
问题 ; 李 伟 华 等 基于 非 饱 和 多 孔 介质 波动 方程 , 研 
究 了 平面 P 波 入 射 下 场地 土 中 的 波 场 问题 。 为 了 更 
加 合理 准确 地 刻画 场地 土 的 黏 弹性 特性 , 刘 林 超 


等 "0 借助 分 数 阶 黏 弹性 本 构 关系 ,在 一 维 波 动 模型 
的 基础 上 研究 了 坚 直 向 上 传播 的 剪 切 地 震波 作用 下 
分 数 阶 茜 弹性 场地 土地 震 放 大 效应 。 段 玮 玮 等 
基于 Biot 动力 固 结 方程 ,将 土 骨架 视 为 黏 弹性 体 ,并 
采用 分 数 阶 导数 本 构 关系 来 描述 土 骨架 的 应 力 应 变 
本 构 关系 ,在 频率 域内 研究 了 饱和 分 数 阶 黏 弹性 十 
层 竖 向 振动 放大 效应 。 很 明显 ,这 些 研 究 是 将 场地 
土 视 为 各 向 均匀 同性 介质 ,对 于 场地 土 内 任意 点 的 
应 力 没有 考虑 周围 各 点 的 影响 ,也 没有 考虑 土 体 颗 
粒 和 孔隙 尺度 等 非 局 部 效应 的 影响 02231 。 
等 5 指出 连续 体内 某 一 点 的 应 力 ,不 仅 只 与 该 点 
的 应 变 有 关 , 还 与 连续 体内 所 有 点 的 应 变 及 变形 历 
史 有 关 , 为 了 能 够 有 效 反映 尺度 因素 的 影响 ,Eringen 
提出 了 非 局 部 理论 。 当 前 非 局 部 理论 主要 被 广泛 应 
用 在 纳米 领域 ,而 在 岩 土工 程 或 其 他 领域 的 应 用 还 
未 得 到 重视 。Li 等 ' 1 首次 将 非 局 部 弹性 理论 和 Biot 
饱和 土 理论 相 结 合 , 提 出 了 一 种 非 局 部 Biot 饱和 土 
理论 ,并 研究 了 流体 饱和 多 孔 介 质 中 波 的 传播 特性 。 
徐 长 节 等 291 在 非 局 部 Biot 理论 的 基础 上 研究 了 饮 
和 土 中 深 埋 圆 柱 形 衬砌 对 平面 弹性 波 的 散射 问题 。 
本 研究 考虑 土 体 非 局 部 效应 的 影响 ,基于 非 局 部 弹 
性 理论 和 分 数 阶 黏 弹性 本 构 理论 ,研究 基于 非 局 部 
理论 的 分 数 阶 黏 弹性 场地 土 振动 放大 效应 ,分 析 非 
局 部 效应 等 对 系统 动力 特性 的 影响 。 
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1 基于 非 局 部 的 分 数 阶 医 弹性 场地 土 
振动 控制 方程 


图 1 所 示 厚 度 为 h 的 场地 土 层 , 场 地 土 下 面 为 
无 穷 大 的 基 岩 ,场地 土 受 向 上 传播 的 剪 切 地 震波 的 
作用 ,地 震波 满足 

v(t) =Doew (1) 
其 中 :wo 为 剪 切 波 横 向 位 移 幅 值 ;w 为 圆 频率 ;i 为 虚 
数 单位 。 


多 1 黏 弹 性 场地 土 


Fig.1 Viscoelastic field soil 


NI 为 了 考虑 士 体 中 任意 点 周围 整个 区 域内 其 他 点 
玉 点 外 应力 的 及 响 运用 Eringen 提出 的 非 局 部 理 
论 来 描述 土 体 的 力学 特性 ,假定 场地 土 为 均 质 线 弹 
向 同 性 介质 ,在 动力 情况 下 ,不 计 体力 时 的 基本 方 
程 为 " 


Cfos = pu 
- ， Cj) |x | = ,e004) CimendV( 广 ) 


1 
Ey = FUiy + i) 


(2) 
其 中 :o; 为 非 局 部 应 力 ;ej; 和 ei 为 应 变 张 量 ;ui 为 位 


移 张 量 ; Cw 为 四 阶 弹性 张 量 ; 积分 |X( 1r-r|， 
ea) 是 非 局 部 模 量 ,eoa 为 表征 非 局 部 效应 的 非 局 部 
核 函 数 ,eo 和 a 分 别 为 非 局 部 材料 常数 和 内 特征 长 
度 ,可 根据 原子 晶 格 动态 波 的 弥散 曲线 拟 合 等 途径 
得 到 , |r -| 为 欧式 距离 。 方 程 (2) 第 2 式 的 微分 
形式 为 


[1- (eva) Vio=0 (3) 
其 中 :cr 为 非 局 部 应 力 张 量 ;o" 为 经 典 应 力 张 量 ,对 
于 场地 土 , 主要 指 土 体 颗 粒 尺 度 等 ;V 为 Laplace 
算 子 。 


由 于 所 研究 的 问题 为 一 维 问题 ,V? = a 
简化 为 
1d _ 下 
| 1 ~ (eo0) Jr = (4) 


其 中 :o, 为 场地 土 非 局 部 剪 切 应 力 ;cr: 为 场地 土 经 
典 剪 切 应 力 。 设 场地 土 为 秋 弹 性 介 且 应 力 应 变 
关系 满足 分 数 阶 Kelvin 本 构 关系 ， 
or=CG | (5) 
其 中 :7 为 场地 土 的 黏 弹 比 ; G0 为 剪 切 弹性 模 量 
y.() 为 切 应 变 ;D" 为 a(0 <a <1) 阶 Riemann-Li- 
ouville 分 数 阶 微分 算 子 ， ee 


a 

D°[x(i)] = i 人 dl 
其 中 TT(，… ) 为 伽 玛 函数 。 

场地 土 受 竖 直 向 上 剪 切 地 震波 作用 时 只 产生 横 
向 位 移 w(z,?) , 则 式 (2) 的 第 1 式 为 

do,(z 2 z， 
其 中 pp 为 场地 土 密度 。 场 地 土 的 应 变 -位 移 为 
y, (三 于 t) (8) 


对 式 (7) 两 端 进行 [1 - (ea)” 45] 运 算 , 并 将 


SN sa , 
Gl 0 ape d 9] 
pl1 -C000) 3] sD -0 (9) 
当 e6=0,a=1 时 , 式 (9) 退 化 为 不 考虑 局 部 效 


应 的 经 典 备 弹 性 解 , 当 e。 =0,a=0 时 , 式 (9) 退 化 为 
不 考虑 局 部 效应 的 弹性 解 。 


2 基于 非 局 部 理论 的 分 数 阶 黏 弹性 场 
地 土地 震 放 大 系数 


由 式 (1) 知 ,地 震波 为 简 谐 地 震波 , 则 上 覆 场 地 
土 位 移 满 足 


u(z,t) =u(z)e™ (10) 
其 中 4(z) 为 横向 位 移 幅 值 。 将 式 (10) 代 入 式 (9)， 
并 考虑 分 数 导数 性 质 , 两 端 消去 e”, 则 有 


CN +7" (iw)" 4 +pw’ [i (eo | =0 


a) dz 
(11) 


G I 
令 ,c. = or We 
po Cs h 
~ 则 式 (11) 无 量 纲 运算 为 
Tu et (12) 
dz” 
3 
其 中 心 = 一 


a 


a 一。 求解 方程 式 (12) 可 得 
1 +i’k ww -7T 


u(z) =AcosAz +BsinAz (13) 
其 中 ,4、B 为 待定 系数 ,可 由 场地 土 上 下 边界 条 件 确 


定 。 忽 略 高 阶 量 的 影响 ,由 式 (3) 可 得 考虑 非 局 部 
效应 的 应 力 为 
Or = ay 十 全 Voy (14) 
由 此 可 得 简 谐 竖 向 地 震波 作用 下 的 前 切 应 力 的 
表达 式 为 


GZ) =A(L+ick" ao")( -4sinAz+BcosAz) (15) 
其 出 ,a.(z) = 5.(z)]Go,6.(z) 为 非 局 部 剪 切 应 力 
网 人。 

刁 如 图 1 所 示 的 弹性 场地 顶 面 的 位 移 和 底面 的 前 
[人 
人 u(z) = z=1 
I z=0 0 


由 式 (13) 式 (15) 和 式 (16) 可 得 


hh? 


VoCOSAZ 


u(z)= (17) 


SZ cosA 
@( ) = ea) Ch nm) 
定义 场地 土 任意 深度 处 的 位 移 幅 值 与 基 岩 位 移 
眉 利 的 比值 为 位 移 放 大 系数 , 剪 切 应 力 幅 值 与 基 岩 
位 移 幅 值 的 比值 为 应 力 放 大 系数 , 则 有 
u(z) 
B= 


Vo 


(18) 


_ lcosAz 
cosA 


(19) 


a sinA z 
) cosA 


|- TA ) (1 +i*k 


(20) 

当 7=0 时 ,问题 的 解 退化 为 不 考虑 局 部 效应 的 

分 数 导数 忒 弹性 解 的 情况 , 当 e。=0,a=1 时 则 退化 

为 不 考虑 局 部 效应 的 经 典 医 弹 性 解 的 情况 , 当 eo = 

0,a=0 时 则 退化 为 不 考虑 局 部 效应 的 弹性 解 的 
情况 。 


3 数值 算 例 分 析 


为 了 研究 场地 土 非 局 部 效应 和 本 构 模 型 参数 等 


对 场地 土地 震 放大 效应 的 影响 ,这 里 以 数值 算 例 的 
形式 对 位 移 放 大 系数 B6 和 z/h =0.5 处 应 力 放大 系 
数 进行 分 析 讨 论 。 未 作 说 明 时 ,参数 的 取 值 分 别 为 : 
Qa =0.5,k=0.8,7=0.05。 图 2 给 出 了 文献 [10] 和 
本 研究 (T=0) 位 移 放 大 系数 的 对 比 , 当 不 考虑 非 局 
部 效应 时 ,本 研究 的 解 与 文献 [ 10] 得 到 的 解 非常 接 
近 , 可 见 本 研究 方法 的 正确 性 。 图 3 ~ 图 8 为 位 移 放 
大 系数 B 和 z/h =0.5 处 应 力 放 大 系数 随 无 量 纲 
频率 hw/c, 的 变化 曲线 ,可 以 看 出 ,在 低频 时 ,位 移 放 
大 系数 和 应 力 放大 系数 随 频率 变化 曲线 存在 波动 。 
图 3 和 图 4 为 非 局 部 效应 对 位 移 放 大 系数 B 和 zy 
=0.5 处 应 力 放大 系数 的 影响 , 当 7=0 时 为 不 考 
虑 土 体 非 局 部 效应 的 情况 ,很 明显 非 局 部 效应 的 情 
况 可 以 退化 到 不 考虑 非 局 部 效应 的 情况 。 在 低频 时 
非 局 部 效应 系数 + 对 位 移 放 大 系数 几乎 没有 影响 ， 
而 在 高 频 时 ,7 越 大 ,位 移 放 大 系数 越 小 ; 而 非 局 部 
效应 系数 对 应 力 放大 系数 & 的 影响 较 大 ,在 低频 时 
应 力 放大 系数 & 随 频 率 变化 曲线 波动 较 大 , 随 着 频率 
的 增 大 逐渐 趋 于 稳定 , 且 7 越 大 ,应 力 放大 系数 越 小 。 


文献 [10] 的 解 
本 研究 解 


0 5 10 15 20 


图 2 文献 [10] 和 本 研究 位 移 放 大 系数 对 比 


Fig.2 Comparison of displacement amplification coefficient 


between Ref. [10] and this paper 
分 数 导数 的 阶 数 对 位 移 放大 系数 B6 和 应 力 放 大 
系数 上 的 影响 较 大 (图 5 和 图 6) , 且 分 数 阶 黏 弹性 
模型 的 解 可 以 退化 到 经 典 舌 弹 性 模型 的 情况 。 分 数 
阶 黏 弹性 场地 土 位 移 放大 系数 B 随 频 率 变 化 曲线 较 
经 典 黏 弹性 模型 波动 要 大 , 且 随 着 分 数 导 数 阶 数 的 
减 小 波动 越 大 ,位 移 放大 系数 B 越 大 ; 随 着 分 数 导数 
阶 数 的 减 小 应 力 放 大 系数 志波 动 更 厉害 。 分 数 导 数 
的 阶 数 是 由 场地 土 的 力学 性 质 决 定 的 , 可见, 场地 土 
的 力学 性 质 对 场地 土 的 地 震 效 应 的 影响 较 大 , 且 分 
数 阶 黏 弹性 模型 可 以 通过 数据 拟 合 土 体 蠕 变 曲线 来 

刻画 不 同性 质 的 场地 土 ,应 用 范围 更 广 。 
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3 非 局 部 效应 对 位 移 放大 系数 的 影响 


Fig.3 Influence of non local effect on the 
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displacement amplification factor 
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加 加 无 量 纲 频率 hoic, 
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a== 图 4 非 局 部 效应 对 应 力 放 大 系数 的 影响 
>< Fig.4 Influence of non local effect on the 
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可 stress amplification factor 
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图 5 分 数 导 数 的 阶 数 对 位 移 放 大 系数 的 影响 
Fig.5 Influence of the order of the fractional derivative 
on the displacement amplification factor 
土 体 秋 性 « 的 有 影响 见 图 7、 图 8, 对 于 场地 土 下 
为 基 岩 的 情况 ,上 禾 场 地 土 的 厅 性 对 位 移 放 大 系数 
的 影响 与 频率 有 关 ,频率 较 小 时 , 土 体 黏 性 k 越 大 ， 
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位 移 放 大 系数 越 小 ,而 在 高 频 时 则 相反 , 土 体 和 针 性 < 
越 大 ,位 移 放 大 系数 越 大 (图 7); 土 体 条 性 < 越 大 ， 
应 力 放大 系数 越 大 ( 图 8)。 土 体 深度 zx 对 应 力 放 
大 系数 的 影响 见 图 9, 由 于 土 体 深 度 z/h 越 深 , 越 接 
近 基 岩 , 所 以 土 体 深度 z/h 越 深 处 的 应 力 放 大 系数 
越 大 , 且 深 度 z/h 对 应 力 放 大 系数 的 影响 越 大 。 
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图 6 分 数 导 数 的 阶 数 对 应 力 放大 系数 的 影响 
Fig.6 Influence of the order of the fractional derivative 


on the stress amplification factor 
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7 土 体 条 性 对 位 移 放 大 系数 的 影响 


Fig.7 Influence of soil viscosity on the displacement 
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8 土 体 黏 性 对 应 力 放 大 系数 的 影响 


Fig.8 JInfluence of soil viscosity on the stress amplification factor 
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图 9 土 体 深 度 对 应 力 放 大 系数 的 影响 
Fig.9 Influence of soil depth on the stress 
amplification factor 

人 
4C 结 论 


C 
6 未 研究 在 提 局 部 理论 攻 弹 作 理 论 和 分 数 导 数 
于 演 的 基础 上 对 荔 强 性 场地 二 的 竖 向 振动 放大 效应 
进行 了 分 析 ,并 分 析 了 土 体 非 局 部 效应 力学 性 质 等 
对 地 震 放大 系数 的 影响 。 主 要 结论 如 下 。 

加 1) 低频 时 ,位 移 放大 系数 和 应 力 放 大 系数 随 频 
才 化 发 生 波动, 计时 运 浙 光 于 稳定 ， 

2) 非 局 部 效应 的 情况 可 以 退化 到 不 考虑 非 局 
钥 应 的 情况 ,分 数 阶 攻 弹性 解 可 以 退化 到 经 典 黏 


届 多 。 
><3) 土 体 孔 际 尺度 等 非 局 部 效应 对 场地 土地 震 
次 江 的 影响 不 应 被 忽略 。 
:三 4) 场 地 土 的 力学 性 质 
响 较 大 。 

5) 非 局 部 参数 和 上 覆 场 地 土 的 理性 对 位 移 放大 
系数 的 影响 与 频率 有 关 。 

可 见 , 对 于 场地 土地 震 放 大 效应 的 研究 ,需要 综 
合 考虑 土 体 非 局 部 效应 力学 性 质 、 黏 性 等 的 影响 ， 
以 便 更 加 准确 地 为 抗震 设计 提供 依据 。 
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